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Abstract 

We review some properties of the Einstein "Gauss-Bonnet" equations for gravity — also called 
the Einstein-Lanczos equations in five and six dimensions, and the Lovelock or Euler equa- 
tions in higher dimensions. We illustrate, by means of simple Kaluza-Klein and brane cos- 
mological models, some consequences of the quasi-linearity of these equations on the Cauchy 
problem (a point first studied by Yvonne Choquet-Bruhat), as well as on "junction condi- 
tions" . 
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I. Introduction 



Vermeil [JJ as early as 1917, then Weyl j2] and Cartan [3J, showed that the Einstein 
tensor was, in any dimension, the only symmetric and conserved (that is divergence free) 
tensor depending only on the metric, its first and second derivatives, the dependence in the 
latter being linear (see below and e.g. 4 J for a definition of (quasi)-linearity). 

Lanczos [H] in 1932 found a generalization of Hilbert's Lagrangian which is quadratic 
in the Riemann tensor and yields, via a Euler variation with respect to the metric, a ten- 
sor (the Lanczos tensor) which is, like Einstein's, symmetric, conserved and second order 
in the metric but, in contrast to Einstein's, only quasi-linear in these second derivatives. 
Another important property is that, just like Hilbert's Lagrangian is a pure divergence in 
two dimensions and Einstein's tensor identically zero in one and two dimensions, we have 
that in four dimensions the Lanczos Lagrangian is a pure divergence (a Euler density, or 
"Gauss-Bonnet", term 0] [S|) and that the Lanczos tensor is identically zero in four or less 
dimensions (a property already known to Bach 

Lovelock [H] generalized these results in 1971 and obtained, for any dimension, a formal 
expression for the most general, symmetric and conserved tensor which is quasi-linear in 
the second derivatives of the metric and does not contain any higher derivatives. He also 
found the Lagrangian from which that tensor derives by a Euler variation with respect to 
the metric : in dimension D it is, in fact, a linear combination of the [D/2] dimensionally 
continued Euler densities ; hence, in dimensions 5 and 6, the explicit form of the Lovelock 
Lagrangian reduces to a linear combination of Hilbert's and Lanczos' Lagrangians ; in di- 
mensions 7 and 8 its explicit expression as well as that of the derived tensor (which are cubic 
in the curvature) were given by Miiller-Hoissen jHj ; an expression for the developped form 
of the quartic Lagrangian, which is the highest to come into play in 9 and 10 dimensions 
can be found in ref PU] . Given its properties, the Lovelock Lagrangian is therefore the most 
natural generalization of Hilbert's to describe pure gravity in higher dimensional spacetimes. 

The idea now that space-time may have more than four dimensions has been recurrent in 
unified field theories since the original proposal by Kaluza [UJ and Klein [T2|, (see e.g. [J3] for 
an account of the "classical period" , and Jl] for a more recent perspective) — and since the 
renewal of string theories (see e.g. [T^]), it is considered as almost a fact. As for modifying, in 
four dimensions, the Hilbert Lagrangian by the inclusion of terms that are non-linear in the 
curvature, it is an idea that goes back to Weyl |2] and Eddington (THj ; of course, from the 
theorems mentionned above, the Euler-derived tensors (symmetric and conserved) contain 
(in order not to be identically zero in four dimensions) terms in derivatives of the metric up 
to the fourth. In the seventies and early eighties such quadratic Lagrangians were exploited 
in view of renormalizing the quantized theory of linearized general relativity (see e.g. \Tj\ 
for a review of that period) as well as to renormalize the stress-energy tensor of quantized 
matter fields in classical, curved, backgrounds, see ^Hj for a review. They however made 
their most forceful entrance when it was shown that they should arise from string theories, 
see [To] . 

The Lovelock tensor differs however from the tensor derived from a generic non-linear 
correction to the Hilbert Lagrangian in that it contains derivatives of the metric of order no 
higher than the second. The main consequence of this property is to avoid the appearance, 
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by increasing the order of differentiation of the field equations, of new classes of dynamical 
solutions which do not approach the unperturbed, Einsteinian, solutions when the non- 
linear terms in the Riemann tensor tend to zero. Hence the Lovelock tensor is of little use 
in achieving goals which ultimately relie on the possibility of singular perturbations such as 
the renormalization of the graviton propagator ^H] or the type of inflation first discussed 
by Starobinski [20J. However the same property guarantees that the quantization of the 
linearized Lovelock theory is free of ghosts and it was argued that, for this reason, the 
Lovelock Lagrangian should appear in the low-energy limit of string theories (see j2Ij, and, 
eg- EU). 

For all these reasons, gravity theories based on the Lovelock tensor were extensively in- 
vestigated in the late 80 's, starting with the work of Madore j2Sj and Miiller-Hoissen [231 
(see e.g. [2E] for a review). In particular, the consequences of the quasi-linearity of these 
equations on the Cauchy problem and structure of characteristics were studied by Choquet- 
Bruhat in j2H| (see also [2Z1)- The closely connected problem of the wave propagation was 
examined in [2Ej> as well as the resulting difficulties in setting up a Hamiltonian formal- 
ism j2HJ- As for the linear stability of various candidate ground states, see e.g. [SHI- Finally, 
the general setting for establishing conservation laws was given in [3*T] . 

Cosmological models also became a prime focus of interest as the early universe appeared 
to be a privileged arena where unified theories could be probed (see e.g. |2B] [23j E2] [33] for 
a taste of the results then obtained, and [33| when gravity is coupled to various other gauge 
matter fields). Interest in the subject then faded away, partly perhaps because priority was 
given to studying the would-be observable properties of the string-theory-predicted dilaton— 
which are better described by four dimensional scalar-tensor theories of gravity (sec e.g. |35j). 

We are however at present witnessing a comeback of pure gravity theories based on 
the Lovelock (or, rather, Lanczos) tensor which is motivated by the invention of "brane 
scenarios" j3E] and, in particular, of cosmological models in which the observable universe is 
described as a four dimensional singular surface, or "brane" , of a five dimensional space-time, 
or "bulk" (see e.g. jSZj and references therein). For an introduction to the current literature 
on Gauss-Bonnet gravity in braneworlds, see [38-55]. 

In this contribution we shall review some properties of the Lovelock field equations — 
more specifically the five dimensional Einstein-Lanczos equations (very often called Einstein 
Gauss-Bonnet) and focus on some consequences of their quasi-linearity, that we shall illus- 
trate by "paradigmatic" Kaluza-Klein or brane cosmological models. We shall also review a 
closely related question which has been debated recently, that is the generalized "junction 
conditions" [5EJ m Lanczos gravity. 



II. The Lanczos and Lovelock Lagrangians and tensors 

a. Quadratic and Lanczos-Gauss-Bonnet Lagrangians 

Consider, in some coordinate system x A , the quadratic Lagrangian density 

V^L q ^^- g (s 2 + Pr LM r LM + 1 R LMNP R LMNP ) (II.l) 

where g is the determinant of the metric coefficients qab, and where R A bcd = ^c^bd ~ ■■■ 
are the components of the Riemann tensor, T^ D being the Christoffel symbols, all indices 
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being raised with the inverse metric g AB ; r B D = R A bad are the Rfcci tensor components, 
s = g BD rBD is the scalar curvature, and (3 and 7 are some constants. Varying this density 
with respect to g AB is a standard calculation which gives 



H AB 5g 



AB 



-gD c V 



c 



{11.2) 



where 



H AB = 2 1 R ALNP R B LNF - 2(2 7 + (3)r^ D R u BAC - ^r u A r CB + 2sr AB - -g AB L q + 



+((3 + Aj)nr AB + -(4 + (3)g AB ns -{2 + (3 + 2j)D AB s 



(II.3) 



and 



V c = \(A 1 + f3)D c r BU + - (4 + f3)g BU D°s - 2D B [((3 + 2 7 )r°" + g uu s \ 5g BD + (II.4) 



+ [2s (g AB g uu - g AU g BU ) + (3(2r AB g uu - r BU g AC - r AC g BU ) + A 1 R BACU }D A 5g BD 

D being the covariant derivative associated with g AB and where □ = g CD DcD- Note that V c 
is defined up to the addition of the divergence-less vector 2c U c with c an arbitrary constant 
and 

tC — ^ n \f m AB„CD BC „AD , AD BC JOD „AB\ 



U u = -D A [(r AB g^ - r B ^g AU + r AU g^ - r^g AB )5g BD 
D A [(r AB g CD -r DC g AB )5g BD ) 



(11.5) 



-(D B a DC ) Sg BD + (r AB g^ - r B( ~' g AB )D A 8g 



AB CD 



BC „AD\ 



BD 



- l>{ - — \sg DC being the contravariant components of the Einstein tensor. 1 



a = r 



The curvature tensor components being second order in the derivatives of the metric 
coefficients, the Lanczos Lagrangian, Lm, is defined by the non trivial combination f3 = —4, 
7 = 1, such that the corresponding Lanczos tensor (of mixed components denoted o~n-\ B ), is 
only second order in the derivatives of the metric coefficients [H] : 

L (2) =s 2 - Ar LM r LM + R LMNP R LMNP (11.6) 

~ \^L m . (II.7) 
Since H AB and crA\ B derive from a Lagrangian they are conserved : 

D A H A = , D A af 2)B = (11.8) 

for any metric g AB (as an explicit calculation, using the Bianchi identities, confirms). As for 
the boundary term (4) it can be written, upon a proper choice of the divergenless vector (5), 
as [S3 

Vg = -AP ABCD D A 8g BD (11.9) 

1 Recall that the variation of the Hilbert Lagrangian density y/— g s with respect to the metric coefficients 
is 5(y/—g s) = y/—g <jab 5g AB + \f—Q DcVfiL where g A b = tab — ^sg A B are the covariant components of 
the Einstein tensor, and where = (g AB g CD — g AC ' g BD )D A 8gBD is uniquely defined. 



„A 
a (2)B 



tjALMN D 

ti- tiBLMN 



2r LM R A 



LBM 



2r AL r BL + sri 
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with 



pABCD _ R ABCD + ^AD g BC _ r BD g AC + r BC g AD _ r AC g B D) _ * {g AB g CD _ g AC g BD\ 



A pedestrian way to see (and define) the quasi- linearity of (which has to be taken 
anyhow when one goes to practical calculations) is to introduce a (local) Gaussian normal 
coordinate system (x A = {w,x^}) such that the line element reads (we use the mostly + 
signature) 

ds 2 = e dw 2 + 7 Mi , dx^dx" . (11.10) 

The surface w = Const, is time- or space-like if e = +1 or —1 and, in terms of the components 
of its extrinsic curvature 

K ~ s -lit < IL11 > 

the components of the Riemann tensor decompose into the gaussian normal version of the 
Gauss-Codazzi-Mainardi equations, that is : 

dKuv 

R-wpwu Q~yj ^P u ^^ ' Rwpvp ^u^pp V * pKjxu 

R\pu P = R\fj,up + ^[K^Kxp — K PP K\ U ] (11.12) 

where V p and R^ vpu are the covariant derivative and the Riemann tensor associated with 
the metric 7^, all Greek indices being raised with the inverse metric 7^. It is then an 
easy exercise to see that, first, u^ B does not contain terms in (dK£ / dw) 2 (as one would 
naively expect since o~^ B contains terms in (Rlmnp) 2 ), second, there are no terms linear in 
(dK^/dw) in and and, third, 

pv _ ivrpau/3 9K a/3 

with 

N paup = 2 {K.K - K 2 )(Y a l vP - Y v l aP ) + A(^K^K v p + ^ v K a p K^)- (11.13) 

-4K{Y u K af3 + ^K^) + SKK^Y 13 ~ %K w K a pl v<i + 4{K^K af3 - K^K V ^) - AeP^ 

where pv v P a i s the restriction onto the surface w = Const, of p ABCD defined in (9); where 
K.K = KpK@ and where the (...) stand for terms of zeroth order in (dK^ /dw). Since 
jypav/3 con tains not only the metric but also its first derivatives with respect to w, the 
Lanczos tensor is, contrarily to Einstein's, only gt/asz-linear in the second derivatives of 
the metric coefficients. 2 We note for further reference that, using the Leibniz rule, the above 
expression can be simplified into 



= 4 ^ - \% J ~ tPU KPa } + - ( IU4 ) 



2 Recall that in Gaussian normal coordinates the components cr™ and a™ of the Einstein tensor do not 
depend on the w— derivatives of the extrinsic curvature and that a^ u — (^ Q <-^ — ^pv^afi^ d ^ fi + ... = 

QK^ r^pv(yK_ _|_ ^ j g n near m them, in the sense that the projector (j^ a ^ y P — y^-y"/ 3 ) contains the metric 

coefficients only 0). 
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with 

jv» = --K^KprK™ + h<K^K u p + -K^iK.K - K 2 ) . 

b. The Lovelock generalisation 

In order now to generalize these results and gain further insight into them, one notices 
that the Hilbert and Lanczos Lagrangians can be rewritten as 

W) = s = R Jl \i 2 

= l^AtA R Jl \i. R Jz \u ( IL1 5) 

where Sj^'fj* is the Kronecker symbol of order 2p : Sj 1 ^" 1 ^ = det 5j, I and J standing for 

Ji, ...lip and Ji, ...Jzp (thus, for example : <5j* j 2 = S^SjI — SjfcSfy. As for the Einstein and 
Lanczos tensor, they can also be rewritten in a similar fashion 

A — A 1 ?A _ _^xAhh r>Jih 

a (l)B — r B 2 B ^BJ X J 2 n hh 

n A _ Ahh hh pJiJ 2 phh fTTIfi^ 

a (2)B - g BJ 1 J 2 J 3 .h K hh K hh- 1,11.10 J 

Lovelock jH] hence generalized the above definitions and showed that the Euler variation 
of the Lagrangian density y/—g L( p ) with 

k P ) = \v 5 tX RJl \n - rJ2p ~ iJ2 ' ( IL1? ) 

was, up to a pure divergence term : 

A _ J x Ah...hp DJ1J2 p^2 P -i^2p /'TT1fi\ 

a (p)B - ~7^TT °BJ 1 ...J 2p K hh ■■■ K hp-ihp- U1.18J 

The developped expressions for L( 3 ) and u^ B were obtained in j2] and an explicit form for 
L( 4 ) can be found in [HI] . 

The Kronecker symbol ".J2 * s zero ^ two or more °fj sa Y) its upper / indices are 
the same ; it is therefore identically zero in space-times of dimension D < 2p ; hence we 
recover that the Einstein tensor is identically zero in one or two dimensions, and see that 
the Lanczos tensor (7) vanishes identically in four and less dimensions (a property first 
discovered by Bach |7.). The Lovelock Lagrangian and tensor are then defined, in a D- 
dimensional space-time, as 

L[D}= E apA^Lfr) , af D]B = £ a p \ 2 ^af p)B (11.19) 

0<P<D/2 0<p<D/2 

where we have set L( ) = 1 and u^ B = —t]S b , where A is a length scale (e.g. the Planck 
scale) and a p are dimensionless parameters, which, for lack of a metatheory or observations, 
have for the time being to be left unspecified. 
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c. The Lovelock Lagrangian and tensor in Cartan's formalism 

The Lagrangians and tensors cr^ B can also be expressed in terms of Cartan moving- 
frame formalism otherwise known, in dimension four, as the Vierbein or tetrad formalism, 
see e.g. |3J. The notation is adapted from [31]. As usual we first define an orthonormal 
frame, that is, a set of D one-forms 9 A which are linear combinations of the differential 
forms dx A naturally associated with the coordinates x A , and chosen so that the metric can 
be expressed in the form 

g = g AB dx A ® dx B = rj AB 6 A <g> 6 B (11.20) 

where tjab = ( — 1, +1, +1) and where the product can be thought of as the symmetric 
tensor product. 

In place of the Christoffel symbols, we then introduce a connection, that is D 2 one-forms 
uj A b which define the covariant derivative by 

D6 A = -u A B ® 6 B (11.21) 

which is extended to arbitrary forms by the Leibniz rule. The connection is torsion-free and 
metric, conditions which can be imposed respectively by the equations 

d9 A + u A B A B = 0, n ac u c B + VBC co c a = 0. (11.22) 

The product here is the exterior (anti-symmetrised tensor) product of two forms and d 
denotes the exterior derivative of a form. 

The curvature 2-form (Cartan's second structural equation) 

Vt A B = du A B + u A c A u c B (11.23) 
defines the components of the Riemann tensor by the expansion 

v a b = \R a bcd e c e D (ii.24) 

where the exterior product symbol A is from now on omitted. 

An important operation is the duality map. We define for each integer p the (D — p)-form 

where ei 1 ,,j D is the completely antisymmetric tensor normalized to 6o...d-i = 1- The Hilbert 
action is then an integral over space-time of the D-form 

£ {1) = n AB 6 AB . (11.26) 

The product here is as always the product in the algebra of forms. As discussed above, for 
each p, the Lovelock Z)-forms 

£ W = m^ei... hp (11.27) 
have similar properties. By Q p we mean p (exterior) products of the curvature 2-form. 
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The Lovelock (D-l)-forms E^a are derived by variation of Cm with respect to the frame. 
Using the Bianchi identities (obtained by exterior differentiation of (23)) one sees that this 
amounts to simply removing one of the frame factors. Hence : 

E {p)A = ^y^6\ h ... l2v . (11.28) 

The components of the Einstein tensor oab are then given by 

E m = -\n BC e* ABC = a AB e* B (11.29) 

(where 0* A = ri AB 6* B ). Variation of the frame yields the same result as variation of the 
metric. This is so because we have forced the torsion to vanish. 

In dimension D = 2p, C( p ) is the Euler density whose integral over the (compact) manifold 
is the Euler number [3] p. In dimension D > 2p, C^, C^, sometimes called "dimension- 
ally continued" Euler densities, become "dynamical" in the sense that their Euler variations, 
to wit the Lovelock forms Em A , E/ P ) A are no longer identically zero. 

For each term (27) we introduce the Sparling forms <J( p )a and t^a defined as 

cr {p)A = -^ C ^ 3 .../ 2p (^Y 3 -^ (H.30) 

t (p)a = {u A B uJ CD e BCDh ... l2p - .V% C ; r ,jr 1 )' 3 - l2p . (H.31) 
It can be shown that 

E( p )a = T ( P )A — dcT( p )A ■ (11.32) 

In form formalism this result is due to Sparling. A proof can be found elsewhere [31] . In the 
case p = 1 the 3-form and the 2-form can be expressed in terms of Christoffel symbols and 
become respectively the 'pseudo-tensor' and the 'pseudo-potential'. Equation (32) becomes 
the conservation of the pseudo-tensor in vacuo. 

Once one has written the Lovelock action in frame formalism it is an easy task to show 
that the field equations are of second order but, apart in the D < 4 case, only quasi-linear, 
in that the coefficients can contain derivatives. We have shown in fact in (31) that each 
contribution to the Sparling (D — l)-form t a is the product of a term quadratic in the 
connection and a power of the curvature. Equation (32) states that the field equations are 
equivalent to the condition that this form be closed. Since the exterior derivative of the 
curvature factors do not contribute, because of the Bianchi identities, the equations are of 
second order, except for the limiting value of p, in which case they become identities. Because 
however of the curvature factors the chacteristic surfaces can be more complicated than a 
simple light cone [2H| |2Tj 3 . 

In the following we shall content ourselves with simple examples taken from cosmological 
models showing the kind of problems which may arise because of the quasi-linearity of the 
Lovelock Lagrangian and tensor. 

3 For example, the propagation of high frequency metric perturbations h A B — £AB& lip such that 
D AB h C D ~ £,A£,Bh C D with £ A = d A f and solution of S<j^ p}B = (N.B. : not Y, p a p \ 2( - p ~V 5<j^ b = 0) 
is determined by the vanishing of the determinant : A = C 2 ^(p-i) cr ^, i !i)?A^B. (see ref (8) in [23) 
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III. Examples from Kamza-Klein Cosmology 



Consider the situation when the Z)-dimensional manifold is the product of an "external" 
(d + l)-dimensional Friedmann-Lemaitre-Roberston- Walker spacetime and a n-dimensional 
compact "internal" space of constant curvature. In an appropriate coordinate system the 
line element can be written as 



ds 2 



—dt 2 + g^ v dx^dx v + g a \ ) dx a dx h 

= -dt 2 + a 2 {t)dQ 2 + a 2 n {t)dQ 2 n (III.l) 

where dQ 2 , n describe (d, n)-dimensional maximally symmetric spaces. It is then a matter 
of straightforward calculation to obtain the components of the Lovelock tensor (11.19). We 
shall write them as 



>[D]0 



If 

2 







2 v 



, . a d 

Id H 9d 

ad 



a 
a„ 



a 



[D]b 



~2 6b 



In H h n H g n 

ad a n 



(III.2) 



where the coefficients F, fd, n , 9d,n, hd, n are polynomials in hd, n = ^dd, n /a dyn and Ad yU = 
^ 2 kd,n/ a dn + with k d , n taking the values (+1, 0, —1) and a dot representing d/dt. Their 
explicit expressions have been obtained up to cubic order by Miiller-Hoissen j^]. They can 
be generalized to any order when either both spaces are static or both spatially flat (see 
Deruelle and Farina-Busto in [B2])- The field equations 



0i 



D]B 



T, 



B 



(III.3) 



with Tg the components of the stress-energy tensor of a perfect fluid, read 



O'd , Q>n 7 

— 9d H h c 

ad a n 



ifd + VdF) 



O'd 7 Q"n 
"n ~r g-n 

ad a n 



ifn + »nF) (III.4) 



where T 00 = p is the energy density of matter and fid,n its adiabatic indices in the external/ 
internal spaces (defined by = fi e p5^, T£ = p e p5l). 

These equations can serve as a basis to study various vacuum (p = pd,n = 0) states and 
their stability. For example, if we impose the external space to be Minkowski spacetime (d = 
3, kd = 0, Xd = 0), then the coefficients a p must satisfy two constraints (see e.g. j2H] Pi] [3^] ) 



E 



nl 



a,. 



0<p<D/2 



[11 



2p)\ al p 



E 



a,. 



(n-1)! 



0<p<D/2 



[n 



5l = o 



(III.5) 



The first gives the radius of the internal space in terms of its dimension n, curvature k n 
and the parameters a p ; the second gives a , say, in terms of the other a p . Of course the 
solution may not be unique. The zero mode linear stability of these ground states, that is, 
the behaviour of the solutions of (4) when expanded at linear order, was studied in [32j and 
shown to impose further constraints on the parameters a p . (Linear stability of the other 
modes was examined by Ishikawa jBOj-) 
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The full-fledged equations (4) can also serve to build cosmological models. Indeed, when 
d = 3 and the size of the internal space is small enough, such a geometry is a potentially 
realistic candidate for describing today's universe. However the (numerical) analysis of the 
dynamical evolution from an initial ground state where both the internal and external spaces 
are of Planckian size to a final state where the external universe evolves in a quasi-Friedman 
way and the internal space freezes out, has not been performed. An important word of 
caution is that if the sign of the determinant 



is positive, say, near the initial state and negative near the final state, then no dynamical 
evolution from one to the other is possible since the system (4) exhibits a breakdown of 
predictability when the determinant goes through zero. 

To illustrate this kind of possible pathology due to the quasi-linearity of the equations, 
consider as an example the case d = 4 (with fc 4 = 0), n = 0. Then the field equations reduce 
to the quadrature 



As studied by Deruelle and Farina-Busto universes obeying such an evolution law can 
exhibit pathological behaviours such as "coming into being" at a time t = —t\ which is not a 
curvature singularity and then, the solution of (7) being multivalued, having the possibility 
of either end up at a curvature singularity at some t = —to or into "nothingness" at t = +ii 
without any divergence in the curvature invariants signaling the approach of that event. The 
origin of such a pathology lies in that the numerator of (7), that is, the coefficient of a in the 
field equation which determines the evolution of the scale factor, vanishes at ±ii. Hence, 
starting from some initial data, at t — for example, the equation for a(t) cannot predict 
its evolution beyond t = t\. (See also, e.g. [21], [TO] and, in the more recent context of 
braneworlds, the "Class I" solutions of [4*§].) 

These pathological cosmological models are just particularly simple examples of the non- 
invertibility of the operator [a^j^j^ ~ 7^7 Q/3) + a 2 N^), with N^ auf3 defined in (11.13). 
We note that rewriting the second order terms in the Lanczos tensor under the form (11.14) 
is of not help to solve the "Cauchy problem" since the first derivative terms hidden in the 
(...) do not combine in an expression proportional to the term in {}. 

IV. Examples from "brane worlds" 

A geometer can construct a "braneworld" as follows : consider a 5-dimensional spacetime 
V + with an edge (an easy way to vizualise this is to imagine, say, the outside surface of a 
half 2-sphere) ; make a copy VI of V + and superpose the copy and the original spacetimes 
onto each other along the edge (this is the so-called Z 2 symmetry) ; one thus obtains a 
spacetime, or "bulk", V5, without an edge (e.g. the outside and inside surfaces of a half 
2-sphere) but which possesses a singular surface, or "brane" £4 whose extrinsic curvature 
is discontinuous : the extrinsic curvature of S 4 embedded in V_ (e.g. the interior of the 
2-sphere) is the opposite of the extrinsic curvature of S 4 embedded in V+ (e.g. the exterior 
of the 2-sphere). 



Det = g n g d - h n h d 



(111.6) 




(111.7) 
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Suppose now that the curvature of V + satisfies the vacuum Lanczos- Gauss-Bonnet equa- 
tions everywhere, except on the edge, that is, is such that 

af 2]B = A5^ + a^ + a af 2)B = (IV. 1) 

the Einstein and Lanczos tensors being defined in (II. 7) (11.16). 

Suppose also that V+ is an anti-de Sitter spacetime everywhere but on the edge. Then, 
because of maximal symmetry, 

Rabcd = —-^{9ac9bd — QadQbc) (IV. 2) 

with 



1 1 / 4aA\ 

T^Ta ( l± f + — J < IV ' 3 > 

because of (1). (One usually chooses the — sign so that, when a — > 0, C? — > — 4, that is the 
Einsteinian value.) 

Finally, suppose, for the time being and the sake of the example, that £4 is flat. To 
find the appropriate edge of V+ (which is locally AdS), one describes it in quasi-conformal 
coordinates such that its line element reads 

ds 2 = dw 2 + e- 2w/ % u dx"dx u (IV .4) 

(C > 0) and one keeps only the w > part, in order that the edge w = be flat. One then 
describes V- with w < 0, so that its line element reads ds 2 = dw 2 + e +2w ^ c r]^ l/ dx tJ 'dx v . 

Consider now the complete space-time V5 U £4. It is locally anti-de Sitter and solution 
of (1) everywhere but on £4. In the chosen coordinates its line element is 

ds 2 = dw 2 + e- 2lwl/ % u dx"dx u . (IV.5) 



The extrinsic curvature of £ 4 is discontinuous (it is K^ v = hj]^ on the V + side and K s 
~\ r \ii,v on the V_ side), that is 

K%=2%S{w) (IV.6) 

where the sign distribution S(w) is +1 if w > 0, and -1 if w < 0. Some components of 
the Riemann tensor therefore exhibit a delta-type singularity (since S'(w) = 25 (w)) and one 
expects that V5 U £4 satisfies the Lanczos-Gauss-Bonnet equations everywhere — that is, £4 
included — but in the presence of "matter" localised on £4, i.e. that one has 

4 ]b = t£v(w) (IV.7) 

where Tg is interpreted as the "stress-energy" tensor of matter in the brane and where T> 
is a distribution localized on £4, i.e. proportional to some linear combination of the Dirac 
delta distribution and its derivatives (see e.g. j3j). 

The question now is to express this "stress-energy" tensor in terms of the discontinuity 
of the extrinsic curvature. 
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Since the line element (5) is written in Gaussian-normal coordinates, the (4 + 1) decom- 
position (11.12) of the components of the Riemann tensor, together with the expression (6) 
for the extrinsic curvature give, near E4 



R 



WflWV 



25(w) + S 2 (w) 



S 2 {w) 
C 2 



{VfipVva-'nuaVvp)- (IV. 



Note that S 2 (w) = 1 everywhere but on w = and is not straightforwardly defined in a 
distributional sense (it is however not "dangerous" when multiplied by just a constant). As 
for the term proportional to S(w), it is distributionally very defined. The Einstein tensor 
reads 



(J. 



6- 



,S 2 (w) 



C 2 



a. 







-6 5£ 



5(w) S 2 (w) 



C 



C 2 



(IV.9) 



Hence, in pure Einstein theory (a = 0, — —A) one recovers the well-known result [HE] 



T. 



W rpyj 



= 5* 1 



(IV.10) 

which is nothing but the "junction conditions" [SH] applied to the problem at hand. 

Now, the Lanczos tensor is only quasi-linear in the second derivatives of the metric 
that is contains, when calculated in gaussian-normal coordinates, terms of the K . K . |^ 
type, see (11.13). More precisely, for an anti-de Sitter bulk and flat brane, (11.13) gives 

24 



a 



£3 



5^S 2 (w)S(w) + ... 



(IV.ll) 



where the dots symbolize non dangerous terms. 



Since the product of the Dirac and sign distributions is not straightforwardly defined, 
various proposals have been put forward to give a meaning to (11) [38-55]. 

Some, starting from (11), have defined S 2 = 1 everywhere, including w = (despite the 
fact that, then, S' = and not 28) and hence rewritten it as 



24 

This way, they obtained that the stress-energy tensor of the brane is 

C v V C 2 ) 



(IV.12) 



(IV. 13) 



Others started from S' = 28 and the Leibniz rule S 2 S' ~ |( i ^ 3 ) / (which is what the rewri- 
ting of (11.13) as (11.14) reduces to in this case 4 ). They then followed the same procedure as 
before to write US 3 )' = US 2 S)' = \S' so that (11) becomes 5 



a 



(2)1/ 



(IV. 14) 



4 However, according to some authors, the Leibniz rule may not apply to products of distributions. 
5 Indeed, if </)(w) is a test function, then dw <j}(w){S 3 )' = 20(0) and hence (S 3 )' = 26 (= S'). 
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and hence obtained the following brane stress-energy tensor 

■«5* (l-^ 
C v \ 3£ 

Finally Deruelle and Dolezel [3H] claimed that the fact that (11) is ill-defined signals that 
one may not treat £4 as a "thin shell" and that the "junction conditions" may depend on 
the microphysics, that is the way one represents it as a "thick" shell. Then the stress energy 
tensor of a flat brane was advocated to be 

Tt = ~6$(l-A%) (IV.16) 



£ v \ C 2 , 

A being a constant which encapsulates the microphysics of the brane. 

When the brane is flat, the difference between the various approaches is immaterial 
as it amounts to a renormalisation of the brane tension. But it matters when one treats 
cosmological models. 

To find a Friedmann-Lemaitre- Robertson- Walker surface £4 in a five dimensional anti-de 
Sitter spacetime, one can introduce various coordinate systems adapted to the problem and 
write the line element (1) under, e.g., its static form 

ds> = -(k + r 2 /C 2 ) dt 2 + 1 -^ 2 + r 2 dQ 2 3) (IV.17) 

with k = or ±1 characterizing the curvature of the maximally symmetric 3-space of line 
element dn 2 {3) = d X 2 + fi(d6 2 + sin 2 9 dip 2 ) (/„ = x, h = sin X, f-i = smh X ). The FLRW 
surface £4 is then defined by 



Ifr? _|_ k/a 2 + 1/ C 2 

r = o(T) , t = t(r) with t = a v - (IV.18) 

a(r) being the scale factor and h = a/a. In the Gaussian normal coordinate system first 
introduced by Binetruy et al. [3*7] 

ds 2 = dy 2 -n 2 (T,y)dr 2 + S 2 (r,y)dn 2 (3h (IV. 19) 

where the expressions for n(r,y) and S(r,y) can be found in [2Z|, the equation for £4 is 
simply y — 0. 

The relevant component of the extrinsic curvature of £4 in V+ is computed to be 



K = K* = Kt = K* = -^ = Jh*+± + ± (IV.20) 



a prime representing Introducing the distribution 

1 1 



C 2 {y) C 2 



+ JC 2 S 2 (y)-1 (IV.21) 
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(which is equal to the constant 1/C 2 everywhere but at y — 0) it is straightforward to obtain, 
using either (11.13) or (11.14), the Einstein-Lanczos tensor, which is proportional to the Dirac 
distribution : 

4a 



a 



[2]r 



-6/C%) 



or 



<7\ 



21 T 



-3/C 5(2/) 



1 



4a . 8 «ir2 C 2/ \ 



(IV.22) 



(IV.23) 



Just as in the case of a flat brane seen above, neither expression (22) nor (23) is straight- 
forwadly defined in a distributional sense. If now, as before, one defines S 2 = 1 everywhere 
including £ in either (22) or (23) one gets the following, very different, results 



p ~ 6 1 1 



4a 
Z 2 " 



'h 2 



k 1 



(IV.24) 



or 



eA 2 + 4 + 4 



£ 2 



4a 8a /, 9 1 
-^ + ^r\h 2 + — + — 



C 



C 2 



(IV.25) 



The first result (24), based on (11.13) and advocated by e.g [20], implies that modifying 
Einstein's gravity by the inclusion of the Gauss-Bonnet correction amounts to nothing more 
(in this particular case) than a coupling-renormalisation. On the other hand, the second 
result (25) based on (11.14) and advocated by, e.g. [HI yields a very different cosmological 
evolution (see, e.g. Lidsey and Nunez in 



In the general case, the issue is then whether or not the stress-energy tensor of matter in 
the brane, as deduced from (11.14), that is 



TH 



(IV.26) 



(Kjf being the extrinsic curvature of the brane in V+), which reduces to (25), is correct or 
not. 



Recently Davis |Qj and Gravanis and Willison [HE], using a former result by Myers [HZ] 
(see also [SHI), proposed to settle the issue positively by considering more closely the bound- 
ary terms which arise when varying the Lovelock Lagrangian with respect to the metric. As 
we have seen in Section II the variation of the Hilbert-Lanczos action is 

5S = 5 f d 5 x V^gL m = I d 5 x^ga [ X 5g AB + f d A x V^l^c (V£ + aV&) (IV.27) 

JM JM JdM V ' y ' 

with V/S given, up to a divergenceless vector by, e.g. (II. 9), and given in footnote 1 ; 
7^ is the induced metric on dAi and n c the unit vector pointing outside the boundary. 
If, now, one imposes that the field equations be obtained by keeping the metric coefficients 
alone fixed at the boundary (their derivatives remaining free) , an extra term must be added 
to the Lagrangian. Now, the boundary term can be written as 

V^in c (Vg + aVg) = <y(V=7Q) +7=7^^ + a^V p W p (IV.28) 
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where only the vector W p = W p {5^ lJiV , V r7 5 / -f fJ-u ) depends upon the choice of the divergenceless 
vector defining V£L. The variation of the action can then be cast into the form (keeping £7^ 
and V pSj^ fixed at the boundaries of dM) 

8Sb = S[ / d 5 x^/— ~g L\2] — / d A x\/— 7 Q ) = 
\Jm JdM J 

= I d 5 x^ga^ B 6g AB + [ d 4 xV=lB^5^. (IV.29) 

JM JdM 

The Chern-Simons form Q and B^ can be written in the language of forms, and their explicit 
expressions are jSZj, [HI], [HH] 



= - K llxv + 4a - - -J lvv - P^KP* ) . (IV.31) 



Q = 2K + Aa(J - 2a flv K' xu ) (IV.30) 

Since is nothing but (26), one can argue convincingly that indeed, (26) is the correct 
junction condition. More precisely, taking the total action as S = Sb + / 24 d 4 Xy/—^L m with 
L m {lpv) the brane matter Lagrangian whose Euler variation yields the brane stress-energy 
tensor <5(i/— jL m ) = —^—jT^Sj^, we get (26) as the brane gravity equations. 

Acknowledgements : One of us (N.D.) thanks Joseph Katz and Cristiano Germani for 
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for Theoretical Physics where this work was completed, for its hospitality. 
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